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1.  Introduction 


Modern  smart  munitions  rely  upon  eleetronic  assemblies  (1).  These  assemblies  are  often 
fabricated  from  commercial-off-the-shelf  components  that  are  not  specifically  designed  to  endure 
the  harsh  loading  environment  of  a  gun  launch.  One  method  of  qualifying  these  assemblies  and 
components  is  to  build  the  projectile  prototypes  and  conduct  live-fire  tests  until  all  design  issues 
are  resolved.  This  approach  is  cost  prohibitive.  Another  approach  is  to  simulate  the  launch 
environment  using  air  guns.  While  less  costly,  there  are  still  drawbacks,  including  acquiring  test 
units,  simulating  the  combined  loading  environment  of  both  setback  and  spin,  and  obtaining 
loading  profiles  that  match  both  the  gun  launch  peak  load  as  well  as  load  duration.  Finally,  the 
finite-element  (FE)  modeling  approach  can  be  used  to  construct  numerical  models  of  the 
components. 

While  FE  modeling  allows  numerous  anticipated  loading  scenarios  to  be  examined  prior  to 
actually  testing  hardware,  issues  related  to  the  modeling  process  do  arise.  These  issues  include 
model  fidelity,  numerical  accuracy,  material  properties,  interface  and  other  boundary  conditions, 
and  solution  time,  ft  would  seem  that  in  the  design  of  smart  munitions  comprised  of  electronic 
components,  we  have  simply  traded  cost  issues  associated  with  actual  testing  with  numerical  and 
other  issues  associated  with  the  EE  modeling  process.  To  some  extent,  this  view  is  correct, 
which  is  why  any  proposed  design  is  rigorously  tested  in  an  actual  gun  launch  environment  prior 
to  fielding.  However,  many  of  the  modeling  issues  mentioned  have  been  addressed  by  prior 
numerical  studies  such  that,  using  reasonable  engineering  judgment,  designs  can  be  screened, 
thus  reducing  the  required  number  of  field  tests.  One  issue  that  has  not  been  addressed 
satisfactorily,  though,  is  the  time  required  to  complete  a  numerical  analysis  of  a  highly  detailed 
EE  model. 

ft  is  now  possible  to  include  more  details  in  the  EE  analysis  with  the  enhanced  pre/post¬ 
processing  tools  because  of  advances  in  EE  modeling  tools,  analysis  codes,  and  the  increase  in 
hardware  computational  capabilities.  However,  the  inclusion  of  more  details  typically  leads  to 
longer  analysis  times.  In  fact,  there  is  a  continuing  trend  towards  including  as  much  of  the  “real” 
physics  as  possible  in  EE  models.  This  conflicts  with  a  designer’s  need  for  quick  analysis  turn¬ 
around  times  so  that  designs  can  be  evaluated  and  improved. 

Techniques  such  as  Guyan  reduction  have  been  used  for  static  analysis  to  reduce  the 
computational  cost  associated  with  the  analysis  of  large,  complex  models  (2).  These  techniques 
have  also  been  applied  to  transient  analysis  problems  where  inertial  effects  are  important  even 
though  the  applicability  of  these  reduction  techniques  for  this  class  of  problems  is  questionable. 
Direct  solutions,  whether  by  implicit  or  explicit  solvers,  of  the  full  model  are  thus  often 
employed  to  avoid  the  questionable  application  of  these  numerical  approximations  on  the  time 
dependent  solution,  even  though  computational  analysis  time  for  these  models  is  substantial 
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(3,  4).  The  analysis  time  required  for  these  models  inereases  even  more  dramatieally  when  the 
model  length  seales  are  highly  disparate.  For  instanee,  a  typieal  artillery  round  is  on  the  order  of 
one  meter  in  length,  while  the  eapaeitors  used  in  the  electronic  components  are  on  the  order  of 
millimeters  or  less  in  size.  Inclusion  of  this  small  capacitor  in  the  overall  model  is  not  feasible, 
since  inclusion  of  these  components  in  the  model  leads  to  excessive  computational  times, 
although  it  is  the  survivability  of  the  electronic  components  that  are  of  concern. 

A  possible  resolution  to  this  issue  is  to  develop  a  global/local  modeling  approach  where  the 
response  of  the  small  components,  the  local  model,  can  be  captured  without  necessarily 
including  detailed  FE  models  of  these  small  components  in  the  global  FE  model  of  the  projectile. 
This  approach  would  thus  not  incur  the  computational  penalty  due  to  the  small  time  steps 
required  if  the  detailed  EE  local  models  were  included,  as  in  a  baseline  EE  model.  There  are 
several  variations  to  this  globaFlocal  modeling.  One  approach  is  to  replace  small  components 
and  intricate  structures  with  generic  coarsely-meshed  representations  that  have  material 
properties  that  approximate  the  response  of  the  actual  local  structure  and  occupy  the  same 
volume  (i).  The  boundary  conditions  in  these  coarse  representations  may  have  significant 
errors,  as  well  as  other  issues  such  as  appropriate  stiffness  properties,  mass  densities  and 
rotational  inertia  properties.  A  second  approach  is  to  represent  the  local  model  as  a  rigid  body. 
This  approach  for  the  overall  model  can  solve  some  of  the  issues  previously  mentioned,  but 
introduces  concerns  about  the  behavior  of  stress  waves  at  the  interface  between  the  elastic  and 
rigid  parts  of  the  EE  model.  Still  another  approach  is  to  simply  replace  the  small  components 
with  lumped  masses  and  ignore  stiffness  effects  completely.  While  this  approach  doesn’t 
preserve  the  volume  requirements,  it  does  work  well  for  uniaxial  response.  Einally, 
superelement  formulations  and  Guyan  reduction,  or  similar  techniques,  can  be  used  to  convert  a 
local  model  into  an  equivalent  superelement. 

These  techniques  can  be  grouped  into  static  and  dynamic  reduction  techniques  (5).  The  static 
reduction  techniques,  such  as  Guyan  reduction,  omit  the  inertial  effect  and  focus  on  the  system 
stiffness.  Similarly,  the  dynamic  reduction  techniques  primarily  focus  on  the  inertial  aspect  of 
the  system,  but  can  also  account  for  both  the  inertial  and  stiffness  effects  for  a  range  of 
frequencies.  Among  these  dynamic  techniques,  the  system  equivalent  reduction  expansion 
process  (SEREP)  is  a  technique  that  can  both  reduce  the  number  of  degrees  of  freedom  (DOE) 
and  maintain  the  correct  dynamic  structural  response  for  a  user-selected  frequency  range. 

SEREP  has  traditionally  been  used  to  aid  in  correlating  modal  test  data  and  modal  numerical 
models  for  PE  modeling  validation.  However,  the  approach  appears  to  be  a  viable  method  for 
representing  the  response  of  complex  structures  without  incurring  some  of  the  deficiencies  noted 
for  other  techniques.  Details  of  the  theory  behind  SEREP  can  be  found  in  several  papers  (6).  In 
this  report,  we  describe  a  method  for  implementing  SEREP  into  the  Lawrence  Livermore 
National  Laboratory  (LLNL)  DynaSd  PE  modeling  software  (7)  to  allow  global/local  analysis. 


2 


2.  SEREP  Implementation 


The  idea  of  the  proposed  approach  is  to  replace  portions  of  the  global  structure  with  their  modal 
representations  to  achieve  faster  solution  times  by  effectively  increasing  the  characteristic  time 
step.  The  difference  between  the  SEREP  approach  and  traditional  modal  analysis  is  that,  while 
modal  analysis  solves  the  problem  in  terms  of  generalized  modal  DOE,  SEREP  reformulates  the 
problem  in  terms  of  the  physical  DOF.  Therefore,  a  EE  model  can  be  partitioned  into  sections 
represented  as  traditional  EE  meshes  while  other  sections  are  represented  by  their  SEREP 
approximations.  Since  both  EE  and  SEREP  formulations  are  in  terms  of  the  physical  DOF,  it  is 
straightforward  to  relate  the  global  and  local  DOF.  To  illustrate  and  clarify  this  distinction, 
consider  the  simple  model  shown  in  figure  1(a).  This  model  was  meshed  traditionally,  as  seen  in 
figure  1(a),  and  also  using  a  mixed  traditional- SEREP  representation,  seen  in  figure  1(b). 


Figure  1 .  Baseline  model  with  (a)  standard  FE  mesh  and  (b)  SEREP  reduced  model. 

In  the  standard  modeling  approach,  either  mesh  congruency  can  be  maintained  at  the  interface 
between  the  two  components  such  that  all  DOF  at  the  interface  are  merged,  or  an  interface  can  be 
defined  which  relates  the  DOF  of  one  component  to  the  other,  as  seen  in  figure  1(a).  In  LS- 
Dyna,  or  ELNL  DynaSd,  this  type  of  interface  is  referred  to  as  a  tied  interface.  In  general,  it  is  a 
multipoint  constraint  equation.  Our  implementation  of  SEREP  can  be  viewed,  although  the 
analogy  is  not  exact,  as  an  extension  of  the  multipoint  constraint  formulation  since  the  DOF  of 
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the  global  model  are  related  by  the  SEREP  representation  of  the  local  model  through  the  retained 
mode  shapes.  This  is  illustrated  in  figure  2  where  node  i  and  node  j  are  related  by  the  modal 
properties  chosen  for  the  representation  of  the  local  model.  This  modal  representation  usually 
results  in  a  full  mass  matrix  for  the  generic  case.  The  mass  matrix  is  not  guaranteed  to  be  full  as 
some  combinations  of  retained  DOE  and  modes  may  lead  to  canceling  of  off-diagonal  terms. 


Figure  2.  SEREP  as  a  multipoint  constraint  approach. 

Implementation  of  the  SEREP  representation  of  the  local  model  requires  the  construction  of  the 
global  mass  and  stiffness  matrices.  Prom  Hopkins  and  Minnicino  (5),  the  SEREP  reduced  mass 
matrix  Ms,  where  the  subscript  5  denotes  SEREP  reduced,  is  given  by 

M,^T^M,T  ,  (1) 

where  Mi  is  the  local  mass  matrix  and  the  transformation  matrix  T  is  given  by 

O 
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o 

\_  tm 

where  the  subscript  r  denotes  the  retained  DOF,  m  denotes  the  retained  modes,  and  t  denotes  the 
truncated  DOF.  When  r  =  m,  the  inverse  required  in  equation  2  is  the  usual  matrix  inverse.  For 
the  general  case  when  r^m,  the  inverse  in  equation  2  is  the  generalized  matrix  inverse  (6). 
Substituting  equation  2  into  equation  1  allows  the  global  mass  matrix  to  be  expressed  as 
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where  the  orthogonality  of  the  mass-normalized  eigenveetors  has  been  utilized.  This  resultant 
matrix  is  singular  if  the  number  of  retained  DOF,  r,  is  greater  than  the  number  of  retained  modes, 
m  (6).  The  meaning  of  the  inverse  shown  in  equation  3  depends  on  the  dimensionality  of  . 
The  superseript -T  denotes  the  inverse  of  the  matrix  transpose,  i.e.,  A'^  =  (A^y\ 

This  SEREP  redueed  loeal  mass  matrix  is  used  to  eonstruet  the  global  mass  matrix 

+M,,  (4) 

where  Md  is  the  diagonal  mass  matrix  obtained  from  the  elements  in  the  global  model  (the  d 
subseript  denotes  diagonal).  M  is  subsequently  used  to  solve  for  the  nodal  accelerations  of  the 
global  model 

a^M-^F.  (5) 

M  may  be  singular  because  Ms  is  singular,  in  which  case,  the  normal  inverse  does  not  exist. 

The  solution  of  equation  5  can  still  be  determined  using  singular  value  decomposition  (SVD)  (5), 
but  this  approach  is  not  discussed  in  this  report.  In  addition  to  being  singular  because  of  rank 
deficiency,  the  reduced  mass  matrix  can  also  be  ill-conditioned.  This  can  also  lead  to  numerical 
problems  with  the  SVD  approach  for  solving  for  the  response.  The  results  presented  in  this 
report  are  for  SEREP  models  that  have  been  constructed  such  that  M  is  nonsingular  and  well- 
conditioned.  This  was  accomplished  by  selecting  local  model  DOE  which  were  on  the  boundary 
of  the  local  model  and  had  corresponding  global  model  DOE  associated  with  the  local  model 
DOE. 

Proper  selection  of  the  retained  DOE  depends  upon  the  actual  boundary  conditions  between  the 
global  and  local  models  as  well  as  consideration  of  the  desired  frequency  response  present  in  the 
local  model’s  SEREP  representation.  As  a  rule  of  thumb,  all  local  modes  below  or  near  the 
highest  frequency  of  interest  should  be  selected.  Also,  only  selecting  nodes  on  the  interface 
between  the  global  and  local  models  typically  leads  to  acceptable  results.  Research  is  currently 
being  conducted  to  develop  more  precise  guidelines  based  on  the  modal  characteristics  of  the 
local  and  global  models. 

The  previous  formulation  of  the  reduced  mass  matrix  was  implemented  in  EENE  DynaSd.  A 
simple  flowchart  of  the  implementation  is  shown  in  figure  3. 
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Figure  3.  Simplified  flow  chart  of  SEREP  implementation. 

It  is  seen  that  implementation  in  LLNL  DynaSd  was  fairly  straightforward.  Routines  were  added 
to  alloeate  storage  arrays  for  the  SEREP  variables,  initialize  these  arrays,  and  then  ealeulate 
effective  forces  acting  upon  the  nodes  associated  with  the  SEREP  nodes.  Initialization  of  the 
reduced  mass  matrix,  and  the  computation  of  its  inverse,  are  done  once  during  the  initialization 
phase.  This  information  is  saved  so  that  computation  of  the  nodal  acceleration  during  the 
solution  phase  is  efficient.  The  only  values  that  must  be  recomputed  at  every  time  step  are  the 
effective  internal  nodal-loads  due  to  the  SEREP  effective  stiffness.  These  internal  forces  are 
computed  using 

(6) 

where  A  is  a  diagonal  matrix  of  the  system’s  eigenvalues  andX,.  „  are  the  nodal  displacements  at 
time  t„  of  the  retained  nodal  DOE. 

The  time  step  cycle  (P)  that  the  solver  uses  to  advance  from  t„  to  t„+i  is  shown  in  figure  4, 
together  with  the  modifications  to  the  DynaSd  logic  required  to  include  SEREP  loading. 

A  time  step  in  DynaSd  starts  with  the  determination  of  the  loads  acting  on  the  nodes,  hirst,  we 
discuss  the  normal  time  step  cycle  for  DynaSd.  This  cycle  will  be  assumed  to  start  at  time  to  with 
all  nodal  displacements,  velocities,  accelerations,  and  applied  forces  known.  The  nodal 
accelerations  are  very  simple  to  compute  because  the  Dyna3d  formulation  results  in  a  diagonal 
mass  matrix.  Thus,  the  nodal  accelerations  are  given  by 

a„,,=M-^F„.  (7) 

The  nodal  velocities  and  displacements  are  subsequently  computed  using  a  central  difference 
update  scheme.  Technically,  DynaSd  uses  a  strain  and  stress  rate  formulation.  However,  for  the 
sake  of  brevity,  the  cycle  simply  indicates  that  the  strains  are  computed  after  the  displacement. 
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Then,  depending  on  the  constitutive  model,  the  element  stresses  are  computed.  Finally,  the 
gradient  of  the  stress  provides  the  internal  forces.  This  internal  force  is  subtracted  from  the 
external  applied  forces  on  a  nodal  basis.  Finally,  the  cycle  completes  with  the  accelerations 
again  being  computed. 

Implementing  SEREP  involves  introducing  two  changes  into  this  time  step  cycle.  The  first 
change,  labeled  “1”  in  figure  4,  requires  updating  the  nodal  internal  forces  using  equation  6. 
Entries  in  the  force  vector  that  are  associated  with  SEREP  DOE  are  temporarily  saved  for  later 
use.  DynaSd  then  computes  the  nodal  accelerations  via  the  normal  update  method,  equation  7. 
However,  this  computation  assumes  that  the  mass  matrix  is  diagonal,  which  for  SEREP-related 
DOEs  is  not  necessarily  true.  Consequently,  as  indicated  at  “3”  of  figure  4,  the  actual  nodal 
accelerations  are  computed  using  the  previously  stored  force  vector  and  the  inverse  of  the  mass 
submatrix  that  is  associated  with  SEREP  DOEs.  These  accelerations  are  then  used  to  replace  the 
erroneous  values,  and  the  normal  time  step  cycle  continues.  This  implementation  thus  requires 
minimal  changes  to  the  DynaSd  time  step  cycle,  as  seen  in  figure  4.  The  most  difficult  aspect  of 
the  implementation  was,  in  fact,  determining  how  to  setup  the  appropriate  storage  requirements 
for  the  implementation  to  ensure  that  Dyna3d  storage  requirements  were  not  adversely  affected. 


Figure  4.  LLNL  Dyna3d  time  step  cycle. 

This  implementation  was  initially  tested  using  the  simple  1-D  discrete  spring-mass  model  shown 
in  figure  5.  This  model  was  used  to  validate  the  implementation.  The  local  model  consisted  of 
the  DOE  associated  with  masses  1^,  while  the  global  model  consisted  of  the  DOE  associated 
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with  masses  1  and  4.  Thus,  the  DOF  of  masses  1  and  4  constitute  the  interface  between  the 
global  and  local  models.  The  local  model’s  equations  of  motion  are  given  by 
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where  the  full  mass  of  DOF  2  and  3,  and  half  of  mass  associated  with  DOF  1  and  4  and  the 
springs  connecting  these  DOF,  are  considered  to  comprise  the  local  model.  The  global  model’s 
equations  of  motion,  based  on  the  DOF  associated  with  masses  1  and  4,  together  with  the  springs 
connected  to  the  rigid  walls,  are  given  by 
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(9) 


where  a  step  load,  F,  was  applied  to  mass  1.  The  forces,/;(t)  and  f2{t),  in  equation  9  are  the 
effective  forces  at  the  global/local  interface.  The  eigenvalues  and  eigenvectors  for  the  local 
model  were  computed  using  equation  8  for  free-free  boundary  conditions.  These 
eigenvalues/eigenvector  pairs  were  used  to  construct  the  SEREP  representation  of  the  local 
model.  The  numerical  results  for  the  complete  spring-mass  system  were  computed  as  the 
baseline  response.  Two  SEREP  cases  were  then  analyzed,  hirst,  all  modes  for  the  local  model 
were  retained.  Eor  this  condition,  the  resultant  equations  of  motion  are  identical  to  the  baseline 
system.  In  the  second  SEREP  case,  the  highest  eigenvalue  was  omitted,  so  only  three  modes 
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were  retained  for  the  loeal  model:  the  rigid  body  motion  and  the  first  two  elastic  modes.  Results 
for  the  response  at  masses  2  and  4  for  all  three  cases  are  shown  in  figure  6.  Both  SEREP 
solutions  are  in  excellent  agreement  with  the  baseline  model.  The  effect  of  neglecting  the 
highest  eigenvalue  shows  up  as  a  slight  difference  in  the  response  for  this  case  as  compared  to 
the  baseline  case.  This  simple  spring-mass  model  verified  that  the  implementation  into  DynaSd 
was  correct.  Comparison  of  the  baseline  and  SEREP  results  indicates  that  this  approach 
represents  a  reasonable  method  of  solving  the  globahlocal  modeling  problem.  It  is  important  to 
realize  that  the  current  implementation  of  SEREP  assumes  that  all  modes  below  a  given  upper 
frequency  are  retained.  If  the  system  shown  in  figure  5  was  driven  at  a  frequency  higher  than 
that  of  the  greatest  retained  system  natural  frequency,  then  the  error  between  the  baseline 
solution  and  SEREP  solution  would  increase. 


Figure  6.  Baseline  and  local  results  of  the  discrete  spring-mass  system. 

Next,  the  acceleration  response  of  the  model  shown  in  figure  1(a)  to  a  unit  step  load  was 
computed.  Baseline  results  for  the  full  model  where  all  common  interface  nodes  have  been 
merged  are  shown  in  figure  7.  As  seen  previously  in  figure  1(b),  15  local  model  nodes  at  the 
global/local  interface  are  merged  because  the  meshes  between  the  local  and  global  models  are 
non-congruent.  Therefore,  these  merged  nodes  are  a  subset  of  the  complete  set  of  interface 
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nodes.  The  Y-direction  stress  results  for  the  SEREP  model,  using  these  merged  nodes  at  the 
interface,  are  shown  in  figure  8.  It  is  seen  that  the  results  at  the  interface  are  different,  which  is 
to  be  expected.  However,  the  contour  levels  are  very  similar  a  short  distance  away  from  the 
interface.  Also  note  that  this  test  case  does  not  represent  the  best  correlation  that  can  be  expected 
from  this  technique  since  not  all  interface  nodes  were  used.  Regardless,  the  results  are 
encouraging.  Using  only  the  common  interface  nodes,  the  mass  matrix  that  is  inverted  to  solve 
for  the  nodal  accelerations,  step  “A”  of  figure  4,  is  rendered  nonsingular.  This  is  because  the 
equation  to  be  solved  is  given  by 


F  =  Ma  =  ([Mj  +  [Mj)a, 


(10) 


where  [M^  ]  is  a  diagonal  mass  matrix  obtained  from  the  single-point  integration  formulation  of 
the  EE  elements  at  the  global/local  interface  and  the  \M ^  ]  is  the  potentially  singular  SEREP- 
derived  reduced  mass  matrix.  The  addition  of  these  two  matrices,  for  this  case,  results  in  a 
nonsingular  matrix  that  is  invertible.  This  is  an  empirical  observation,  since  mathematically  the 
resultant  matrix  could  be  singular  for  arbitrary  matrices  \M ^  ]  and  [M^  ] .  If  internal  DOE  were 
selected  from  the  local  model  for  this  case,  then  the  resultant  mass  matrix  is  guaranteed  to  be 
singular. 


Model  fo*"  SCREP 
1  -  5.69995e-03 


Figure  7.  Y  Stress  results  of  3D  model.  Baseline  model  results 
using  only  subset  of  interface  nodes. 
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Figure  3-  Serep  w*lh  only  interface  nodes 
1  -  5.€9986e-03 


Figure  8.  Y  Stress  results  of  3D  model.  SEREP  model  results  using 
subset  of  interfaee  nodes. 

As  was  observed  from  the  Y-direction  stress  global  results,  examining  the  principal  stresses 
shows  that  the  results  are  different  near  the  interface  but  similar  away  from  the  interface,  as  seen 
in  figures  9  and  10. 


Figure  9.  Principal  stress  at  interface — baseline 
model. 
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Figure  10.  Principal  stress  at  interface — SEREP 
model. 

This  observation  allows  for  the  efficient  computation  of  the  high-fidelity  local  solution  by 
applying  the  global  solution  as  local  model  boundary  conditions  where  the  full  and  global 
solutions  appear  equivalent.  In  a  real  structural  analysis,  the  full  system  solution  is  not  known 
and  the  correct  global  boundary  condition  is  an  engineering  judgment  based  on  the  physics  of  the 
problem.  The  appropriate  boundary  conditions  to  apply  to  the  local  model  are  the  time 
dependent  displacements.  Stress/time  histories  could  also  be  used,  but  are  more  likely  to 
introduce  errors  since  stresses  depend  upon  the  displacement  gradients.  While  the  displacements 
of  the  interface  nodes  themselves  could  be  used,  because  of  St.  Venant’s  Principle,  the  localized 
stresses,  strains,  and  thus  displacements  of  these  nodal  displacements  are  not  representative  of 
the  actual  system.  However,  the  use  of  nodes  located  away  from  the  actual  interface  would 
provide  a  better  local  response  where  the  new,  high-fidelity  local  model  now  includes  the 
original  local  model  and  a  part  of  the  original  global  model,  the  size  of  which  is  dictated  by  St. 
Venant’s  Principle. 

For  comparison  with  the  previous  models,  results  for  a  more  representative  model  using  a  tied 
interface  are  presented  in  figure  1 1 .  Unlike  the  SEREP  or  merged  results  presented  previously,  it 
is  seen  that  the  actual  response  at  the  interface  is  smooth.  Again,  this  difference  is  to  be  expected 
for  the  sample  problem  since  not  all  of  the  interface  nodes  were  used  in  either  the  SEREP  or 
baseline  models,  ft  is  hypothesized  that  if  the  local  model  mesh  is  constructed  such  that  all  of 
the  global  model’s  interface  nodes  have  corresponding  nodes  in  the  local  model,  then  a  more 
accurate  response  can  be  obtained  for  both  the  global  and  local  models.  This  is  inferred  because 
the  effect  of  point  loading  on  the  global  model  will  be  reduced  if  the  global  and  local  models  are 
meshed  such  that  the  above  mesh  congruency  is  true. 
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Figure  11.  Effective  stress  at  interface  using  a  tied  interface. 


3.  Future  Work 


Based  on  these  preliminary  results,  the  SEREP  approach  for  global/local  modeling  presents  a 
promising  method  for  reducing  a  EE  model’s  complexity  while  retaining  correct  mass  and 
stiffness  effects.  Remaining  issues  which  must  be  addressed  include  implementing  a  more 
robust  SVD  solver  for  finding  solutions  when  the  mass  matrix  is  singular,  automating  the  task  of 
relating  global  and  local  DOE,  and  automating  the  selection  of  a  proper  set  of  local  DOE  so  as  to 
maintain  accuracy  in  the  global  response  calculation.  None  of  these  issues  are  intractable.  They 
simply  require  diligence  in  formulating  a  reasonable  modeling  methodology  such  that  this 
technique  can  be  a  useful  tool  for  quickly  analyzing  proposed  designs  involving  electronic 
components.  Euture  research  will  examine  SEREP’s  ability  to  use  actual  modal  test  data  to 
generate  the  SEREP  superelement  for  inclusion  in  the  PE  model.  This  would  be  useful  when 
construction  of  a  representative  PE  model  for  a  complex  part  cannot  be  obtained  in  a  timely 
manner,  but  modal  data  either  already  exists  or  can  be  obtained. 
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AMSSB  RCP  SS 

W  NYKVIST 

S  BEAUDOIN 

KANSAS  ST 

NATICK  MA  01760-5019 


1  DIRECTOR 

LOS  ALAMOS  NATL  LAB 
F  L  ADDESSIO  T  3  MS  5000 
PO  BOX  1633 
LOS  ALAMOS  NM  87545 

8  NSWC 

J  FRANCIS  CODE  G30 
D  WILSON  CODE  G32 
R  D  COOPER  CODE  G32 
J  FRAYSSE  CODE  G33 
E  ROWE  CODE  G33 
T  DURAN  CODE  G33 
L  DE  SIMONE  CODE  G33 
R  HUBBARD  CODE  G33 
DAHLGREN  VA  22448 
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1  NSWC 

CARDEROCK  DIVISION 
R  CRANE  CODE  6553 
9500  MACARTHUR  BLVD 
WEST  BETHESDA  MD  20817-5700 

1  AFRL  MESS 
R  THOMSON 

2179  I2TH  STRM  122 
WRIGHT  PATTERSON  AFB  OH 
45433-7718 

2  AFRL 

F  ABRAMS 

J  BROWN 

BLDG  653 

2977  P  ST  STE  6 

WRIGHT  PATTERSON  AFB  OH 

45433-7739 

5  DIRECTOR 
LLNL 

R  CHRISTENSEN 
S  DETERESA 
F  MAGNESS 
M  FINGER  MS  313 
M  MURPHY  L  282 
PO  BOX  808 
LIVERMORE  CA  94550 

I  AFRL  MLS  OL 
L  COULTER 
5851  FAVE 
BLDG  849  RM  ADIA 
HILL  AFB  UT  84056-5713 

1  OSD 

JOINT  CCD  TEST  FORCE 
OSD  JCCD 
R  WILLIAMS 
3909  HALLS  FERRY  RD 
VICKSBURG  MS  29180-6199 

3  DARPA 

M  VANFOSSEN 
S  WAX 

L  CHRISTODOULOU 
3701  N  FAIRFAX  DR 
ARLINGTON  VA  22203-1714 

1  OAK  RIDGE  NATL  LAB 
R  M  DAVIS 
PO  BOX  2008 

OAK  RIDGE  TN  3  783 1  -6 1 95 


I  OAK  RIDGE  NATL  LAB 
C  EBERLE  MS  8048 
PO  BOX  2008 
OAK  RIDGE  TN  3783 1 

3  DIRECTOR 
SANDIA  NATL  LABS 
APPLIED  MECHS  DEPT 
MS  9042 

J  HANDBOOK 
YRKAN 
J  LAUFFER 
PO  BOX  969 

LIVERMORE  CA  94551-0969 

I  OAK  RIDGE  NATL  LAB 
C  D  WARREN  MS  8039 
PO  BOX  2008 
OAK  RIDGE  TN  3783 1 

4  NIST 

M  VANLANDINGHAM  MS  8621 
J  CHIN  MS  8621 
J  MARTIN  MS  8621 
DDUTHINH  MS  8611 
100  BUREAU  DR 
GAITHERSBURG  MD  20899 

I  HYDROGEOLOGIC  INC 
SERDP  ESTOP  SPT  OFC 
S  WALSH 

1155  HERNDON  PKWY  STE  900 
HERNDON  VA  20170 

3  NASA  LANGLEY  RESEARCH  CTR 
AMSRD  ARE  VT 
W  LEBER  MS  266 
F  BARTLETT  JR  MS  266 
G  FARLEY  MS  266 
HAMPTON  VA  23681-0001 

I  FHWA 

E  MUNLEY 

6300  GEORGETOWN  PIKE 
MCLEAN  VA  22 10 1 

I  USDOT  FEDERAL  RAILROAD 
MFATEHRDV3I 
WASHINGTON  DC  20590 
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CYTEC  FIBERITE 
R DUNNE 
D  KOHLI 
R  MAYHEW 
1300  REVOLUTION  ST 
HAVRE  DE  GRACE  MD  21078 

DIRECTOR 

NGIC 

lANG  TMT 

2055  BOULDERS  RD 

CHARLOTTESVILLE  VA 

22911-8318 

3TEX  CORP 
A  BOGDANOVICH 
J  SINGLETARY 
109  MACKENAN  DR 
CARYNC  27511 

DIRECTOR 

DEFENSE  INTLLGNC  AGNCY 
TA5 

K  CREELING 
WASHINGTON  DC  20310 

COMPOSITE  MATERIALS  INC 
D  SHORTT 
19105  63  AVE  NE 
PO  BOX  25 

ARLINGTON  WA  98223 

IPS  GLASS 
L  CARTER 
PO  BOX  260 
SLATER  RD 
SLATER  SC  29683 

COMPOSITE  MATERIALS  INC 
R  HOLLAND 
1 1  JEWEL  CT 
ORINDA  CA  94563 

COMPOSITE  MATERIALS  INC 
C  RILEY 

14530  S  ANSON  AVE 
SANTA  FE  SPRINGS  CA  90670 

SIMULA 
R  HUYETT 
10016  S51ST  ST 
PHOENIX  AZ  85044 


2  PROTECTION  MATERIALS  INC 
M  MILLER 
F  CRILLEY 
14000  NW  58  CT 
MIAMI  LAKES  FL  33014 

2  FOSTER  MILLER 
M  ROYLANCE 
WZUKAS 
195  BEAR  HILL  RD 
WALTHAM  MA  02354-1196 

1  ROM  DEVELOPMENT  CORP 
R  O  MEARA 

136  SWINEBURNE  ROW 
BRICK  MARKET  PLACE 
NEWPORT  RI  02840 

2  TEXTRON  SYSTEMS 
M  TREASURE 

T  FOLTZ 

1449  MIDDLESEX  ST 
LOWELL  MA  01851 

1  O  GARA  HESS  &  EISENHARDT 
M  GILLESPIE 
9113  LESAINT  DR 
FAIRFIELD  OH  45014 

1  MILLIKEN  RESEARCH  CORP 
M  MACLEOD 
PO  BOX  1926 
SPARTANBURG  SC  29303 

1  CONNEAUGHT  INDUSTRIES  INC 
J  SANTOS 
PO  BOX  1425 
COVENTRY  RI  02816 

1  ARMTEC  DEFENSE  PRODUCTS 
SDYER 
85  901  AVE  53 
PO  BOX  848 
COACHELLA  CA  92236 

3  PACIFIC  NORTHWEST  LAB 
M  SMITH 

G  VAN  ARSDALE 
R  SHIPPELL 
PO  BOX  999 
RICHLAND  WA  99352 
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1  ALLIANT  TECHSYSTEMS  INC 
4700  NATHAN  LN  N 
PLYMOUTH  MN  55442-2512 

1  APPLIED  COMPOSITES 
W  GRISCH 

333  NORTH  SIXTH  ST 
ST  CHARLES  IE  60174 

1  CUSTOM  ANALYTICAL 
ENG  SYS  INC 
A  ALEXANDER 
13000  TENSOR  LANE  NE 
FLINTSTONE  MD  21530 

1  AAI CORP 

DR  N  B  MCNELLIS 
PO  BOX  126 

HUNT  VALLEY  MD  21030-0126 

1  OFC  DEPUTY  UNDER  SEC  DEFNS 
J  THOMPSON 

1745  JEFFERSON  DAVIS  HWY 
CRYSTAL  SQ  4  STL  501 
ARLINGTON  VA  22202 

3  ALLIANT  TECHSYSTEMS  INC 
J  CONDON 
ELYNAM 
J  GERHARD 
WVOl  16  STATE  RT  956 
PO  BOX  210 
ROCKET  CENTER  WV 
26726-0210 

1  PROJECTILE  TECHNOLOGY  INC 
515  GILES  ST 

HAVRE  DE  GRACE  MD  21078 

1  PRATT  &  WHITNEY 
C  WATSON 

400  MAIN  ST  MS  114  37 
EAST  HARTFORD  CT  06108 

5  NORTHROP  GRUMMAN 
B  IRWIN 
K  EVANS 
D  EWART 

A  SHREKENHAMER 
J  MCGLYNN 
BLDG  160  DEPT  3700 
1 100  W  HOLLYVALE  ST 
AZUSA  CA  91701 


1  BRIGS  COMPANY 
J  BACKOFEN 
2668  PETERBOROUGH  ST 
HERNDON  VA  22071-2443 

1  ZERNOW  TECHNICAL  SERVICES 
L  ZERNOW 

425  W  BONITA  AVE  STE  208 
SAN  DIMAS  CA  91773 

1  GENERAL  DYNAMICS  OTS 
L  WHITMORE 

10101  NINTH  ST  NORTH 
ST  PETERSBURG  FL  33702 

2  GENERAL  DYNAMICS  OTS 
FLINCHBAUGH  DIV 

K  LINDE 
T  LYNCH 
PO  BOX  127 
RED  LION  PA  17356 

1  GKN  WESTLAND  AEROSPACE 
DOLDS 

450  MURDOCK  AVE 
MERIDEN  CT  06450-8324 

1  BOEING  ROTORCRAFT 
P  HANDEL 

800  B  PUTNAM  BLVD 
WALLINGFORD  PA  19086 

5  SIKORSKY  AIRCRAFT 

G  JACARUSO 
T  CARSTENSAN 
BKAY 

S  GARBO  MS  S330A 
J  ADELMANN 
6900  MAIN  ST 
PO  BOX  9729 

STRATFORD  CT  06497-9729 

1  AEROSPACE  CORP 
G  HAWKINS  M4  945 

2350  E  EL  SEGUNDO  BLVD 
EL  SEGUNDO  CA  90245 

2  CYTEC  FIBERITE 
MEIN 

W  WEB 

1440  N  KRAEMER  BLVD 
ANAHEIM  CA  92806 
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2  UDLP 

G  THOMAS 
M  MACLEAN 
PO  BOX  58123 
SANTA  CLARA  CA  95052 

2  UDLP 

R  BRYNSVOLD 
P  JANKEMS  170 
4800  E  RIVER  RD 
MINNEAPOLIS  MN  55421-1498 

1  NORTHROP  GRUMMAN  CORP 
ELECTRONIC  SENSORS 
&  SYSTEMS  DIV 
E  SCHOCH  MS  V  16 
1745 A  W  NURSERY  RD 
LINTHICUM  MD  21090 

1  GDLS  DIVISION 
D  BARTEL 

PO  BOX  1901 
WARREN  MI  48090 

2  GDLS 
DREES 
M  PASIK 
PO  BOX  2074 
WARREN  MI  48090-2074 

1  GDLS 

MUSKEGON  OPER 
M  SOIMAR 
76  GETTY  ST 
MUSKEGON  MI  49442 

1  GENERAL  DYNAMICS 
AMPHIBIOUS  SYS 
SURVIVABILITY  LEAD 
G  WALKER 
991  ANNAPOLIS  WAY 
WOODBRIDGE  VA  22191 

6  INST  FOR  ADVANCED 
TECH 
HFAIR 
I  MCNAB 
P  SULLIVAN 
S  BLESS 
W  REINECKE 
C  PERSAD 

3925  W  BRAKER  LN  STE  400 
AUSTIN  TX  78759-5316 


1  ARROW  TECH  AS  SOC 

1233  SHELBURNE  RD  STE  D8 
SOUTH  BURLINGTON  VT 
05403-7700 

1  R  EICHELBERGER 
CONSULTANT 
409  W  CATHERINE  ST 
BEL  AIR  MD  21014-3613 

1  SAIC 

G  CHRYSSOMALLIS 

8500  NORMANDALE  LAKE  BLVD 

SUITE  1610 

BLOOMINGTON  MN  55437-3828 

1  UCLA  MANE  DEPT  ENGR IV 
H  T  HAHN 

LOS  ANGELES  CA  90024-1597 

2  UNIV  OF  DAYTON 
RESEARCH  INST 
R  YKIM 
AKROY 

300  COLLEGE  PARK  AVE 
DAYTON  OH  45469-0168 

1  GA  TECH  RESEARCH  INST 

GA  INST  OF  TCHNLGY 
P  FRIEDERICH 
ATLANTA  GA  30392 

1  MICHIGAN  ST  UNIV 
MSM  DEPT 
R  AVERILL 
3515  EB 

EAST  LANSING  MI  48824-1226 

1  PENN  STATE  UNIV 
R  S  ENGEL 

245  HAMMOND  BLDG 
UNIVERSITY  PARK  PA  16801 

1  PENN  STATE  UNIV 
C  BAKIS 

212  EARTH  ENGR 
SCIENCES  BLDG 
UNIVERSITY  PARK  PA  16802 

1  PURDUE  UNIV 

SCHOOL  OF  AERO  &  ASTRO 
CTSUN 

W  LAFAYETTE  IN  47907-1282 
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1  STANFORD  UNIV 

DEPT  OF  AERONAUTICS 
&  AEROBALLISTICS 
STSAI 

DURANT  BLDG 
STANFORD  CA  94305 

1  UNIV  OF  MAINE 

ADV  STR  &  COMP  LAB 
R  LOPEZ  ANIDO 
5793  AEWC  BLDG 
ORONO  ME  04469-5793 

1  JOHNS  HOPKINS  UNIV 
APPLIED  PHYSICS  LAB 
P  WIENHOLD 
11100  JOHNS  HOPKINS  RD 
LAUREL  MD  20723-6099 

1  UNIV  OF  DAYTON 
J  M  WHITNEY 
COLLEGE  PARK  AVE 
DAYTON  OH  45469-0240 

5  UNIV  OF  DELAWARE 

CTR  FOR  COMPOSITE  MTRLS 
J  GILLESPIE 
M  SANT ARE 
S  YARLAGADDA 
S  ADVANI 
D  HEIDER 
201  SPENCER  LAB 
NEWARK  DE  1 97 1 6 

1  DEPT  OF  MTRLS 
SCIENCE  &  ENGRG 
UNIV  OF  ILLINOIS 
AT  URBANA  CHAMPAIGN 
J  ECONOMY 
1304  W  GREEN  ST  115B 
URBANA  IE  61801 

1  UNIV  OF  MARYLAND 

DEPT  OF  AEROSPACE  ENGRG 
A  J  VIZZINI 

COLLEGE  PARK  MD  I9I04 

1  DREXEL  UNIV 
A  S  D  WANG 
3141  CHESTNUT  ST 
PHILADELPHIA  PA  20742 


3  UNIV  OF  TEXAS  AT  AUSTIN 
CTR  FOR  ELECTROMECHANICS 
J  PRICE 
A  WALLS 
J  KITZMILLER 
10 100  BURNET  RD 
AUSTIN  TX  78758-4497 

2  VA  POLYTECHNICAL 
INST  &  STATE  UNIV 
DEPT  OF  ESM 

M  W  HYER 
R  JONES 

BLACKSBURG  VA  24061-0219 

I  SOUTHWEST  RESEARCH  INST 
ENGR  &  MATE  SCIENCES  DIV 
J  RIEGEL 

6220  CULEBRA  RD 
PO  DRAWER  28510 
SAN  ANTONIO  TX  78228-0510 

I  BATELLE  NATICK  OPERS 

B  HALPIN 
313  SPEENST 
NATICK  MA  01760 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRD  ARE  WM  MB 
A  FRYDMAN 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

ABERDEEN  PROVING  GROUND 

I  US  ARMY  ATC 
CSTE  DTC  AT  AC  I 
W  C  FRAZER 
400  COLLERAN  RD 
APG  MD  21005-5059 

87  DIR  USARL 

AMSRD  ARE  Cl 
AMSRD  ARE  O  AP  EG 
M  ADAMSON 
AMSRD  ARE  SL  BA 
AMSRD  ARE  SL  BB 
DBELY 

AMSRD  ARE  WM 
J  SMITH 

AMSRD  ARE  WM  B 
CHIEF 
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AMSRD  ARE  WM  BA 
CHIEF 

AMSRD  ARE  WM  BC 
J  NEWILL 
P  PLOSTINS 
AMSRD  ARE  WM  BD 
P  CONROY 
B  FORCH 
M  LEADORE 
C  LEVERITT 
RLIEB 

R  PESCE-RODRIGUEZ 
BRICE 
A  ZIELINSKI 
AMSRD  ARE  WM  BF 
S  WILKERSON 
AMSRD  ARE  WM  M 
J  MCCAULEY 
S  MCKNIGHT 
AMSRD  ARE  WM  MA 
CHIEF 
L  GHIORSE 
E  WETZEL 
AMSRD  ARE  WM  MB 
J  BENDER 
T  BOGETTI 
J  BROWN 
L  BURTON 
R  CARTER 
KCHO 

W  DE  ROSSET 
G  DEWING 
R  DOWDING 
W  DRYSDALE 
R  EMERSON 
DGRAY 
D  HOPKINS 
RKASTE 
L  KECSKES 
M  MINNICINO 
B  POWERS 
D  SNOHA 
J  SOUTH 
M  STAKER 
JSWAB 
J  TZENG 

AMSRD  ARE  WM  MC 
CHIEF 
R  BOSSOLI 
ECHIN 

S  CORNELISON 
D  GRANVILLE 
BHART 
J  LASALVIA 


J  MONTGOMERY 
F  PIERCE 
ERIGAS 
W  SPURGEON 
AMSRD  ARE  WM  MD 
B  CHEESEMAN 
P  DEHMER 
R  DOOLEY 
G  GAZONAS 
S  GHIORSE 
M  KLUSEWITZ 
WROY 
J  SANDS 
D  SPAGNUOLO 
S  WALSH 
S  WOLF 

AMSRD  ARE  WM  RP 
J  BORNSTEIN 
C  SHOEMAKER 
AMSRD  ARE  WM  T 
B  BURNS 

AMSRD  ARE  WM  TA 
W  BRUCHEY 
M  BURKINS 
W  GILLICH 
B  GOOCH 
T  HAVEL 
C  HOPPLE 
E  HORWATH 
M  NORMANDIA 
J  RUNYEON 
M  ZOLTOSKI 
AMSRD  ARE  WM  TB 
P  BAKER 

AMSRD  ARE  WM  TC 
R  COATES 
AMSRD  ARE  WM  TD 
D  DANDEKAR 
M  RAFTENBERG 
S  SCHOENFELD 
T  WEERASOORIYA 
AMSRD  ARE  WM  TE 
CHIEF 
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1  LTD 

R  MARTIN 
MERE 

TAMWORTH  RD 
HERTFORD  SGI 3  7DG 
UK 

1  SMC  SCOTLAND 
P  WLAY 
DERA  ROSYTH 
ROSYTH  ROYAL  DOCKYARD 
DUNFERMLINE  FIFE  KY  2XR 
UK 

I  CIVIL  AVIATION 
ADMINSTRATION 
T  GOTTESMAN 
PO  BOX  8 

BEN  GURION  INTRNL  AIRPORT 

LOD  70150 

ISRAEL 

1  AEROSPATIALE 
S ANDRE 

ABTE  CCRTE  MDI32 
316  ROUTE  DE  BAYONNE 
TOULOUSE  31060 
FRANCE 

1  DRA  FORT  HALSTEAD 
P  N  JONES 

SEVEN  OAKS  KENT  TN  I47BP 
UK 

1  SWISS  FEDERAL  ARMAMENTS 
WKS 
WLANZ 

ALLMENDSTRASSE  86 
3602  THUN 
SWITZERLAND 

1  DYNAMEC  RESEARCH  LAB 
AKE  PERSSON 
BOX  201 

SE  I5I  23  SODERTALJE 
SWEDEN 

1  ISRAEL  INST  OF  TECHLGY 
S  BODNER 

FACULTY  OF  MECHANICAL 

ENGR 

HAIFA  3200 

ISRAEL 


I  DSTO 

WEAPONS  SYSTEMS  DIVISION 
N  BURMAN  RLLWS 
SALISBURY 

SOUTH  AUSTRALIA  5108 
AUSTRALIA 

I  DBF  RES  ESTABLISHMENT 
VALCARTIER 
A  DUPUIS 

2459  BLVD  PIE  XI  NORTH 
VALCARTIER  QUEBEC 
CANADA 

PO  BOX  8800  COURCELETTE 
GOA  IRO  QUEBEC 
CANADA 

I  ECOLE  POLYTECH 

J  MANSON 
DMX  ETC 

CH  1015  LAUSANNE 
SWITZERLAND 

1  TNO  DEFENSE  RESEARCH 
R  IJSSELSTEIN 
ACCOUNT  DIRECTOR 
R&D  ARMEE 

PO  BOX  6006 
2600  JA  DELFT 
THE  NETHERLANDS 

2  FOA  NATL  DEFENSE  RESEARCH 
ESTAB 

DIR  DEPT  OF  WEAPONS  & 
PROTECTION 
B  JANZON 
R  HOLMLIN 
S  172  90  STOCKHOLM 
SWEDEN 

2  DEFENSE  TECH  &  PROC 
AGENCY  GROUND 
I  CREWTHER 
GENERAL  HERZOG  HAUS 
3602  THUN 
SWITZERLAND 
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1  MINISTRY  OF  DEFENCE 
RAFAEL 

ARMAMENT  DEVELOPMENT 
AUTH 

M  MAYSELESS 
PO  BOX  2250 
HAIFA  3 1 02 1 
ISRAEL 

1  B  HIRSCH 

TACHKEMONY  ST  6 

NETAMUA426II 

ISRAEL 

I  TNO  DEFENSE  RESEARCH 
I H  PASMAN 
POSTBUS  6006 
2600  JA  DELFT 
THE  NETHERLANDS 

1  DEUTSCHE  AEROSPACE  AG 

DYNAMICS  SYSTEMS 
MHELD 
PO  BOX  1340 

D  86523  SCHROBENHAUSEN 
GERMANY 
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